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Abstract 

In 1929 Heisenberg and Pauli laid the foundations of QFT by quantiz- 
ing the fields (method of canonical quantization). This general theory of 
quantized fields has remained undisputed up to now. We show how the un- 
modified Heisenberg-Pauli calculations make sense mathematically by using 
a theory of generalized functions adapted to nonlinear operations. By pro- 
viding an appropriate mathematical setting, nonlinear generalized functions 
open doors for their understanding but there remains presumably very hard 
technical problems, (i) Domains of the interacting field operators: a priori 
the H-P calculations give time dependent dense domains, what is not very 
convenient; (ii) Calculations of the resulting matrix elements of the S op- 
erator: from the unitarity of the S operator as a whole there are no longer 
"infinities," but a priori there is no other hope than heavy computer calcu- 
lations; (iii) Connection with renormalization theory: it should provide an 
approximation when the coupling constant is small. The aim of this paper is 
to present, on the grounds of a standard mathematical model of QFT (a self 
interacting scalar boson field), a basis for improvement without significant 
prerequisites in mathematics and physics. It is an attempt to use nonlinear 
generalized functions in QFT, starting directly from the calculations made 
by physicists, in the same way as they have already been used in classical 
mechanics and general relativity. 
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1 Prerequisites on generalized functions 



As is well known the Heisenberg-Pauli calculations (or ''canonical formalism," 
see [[24l p. 292]), require a concept of generalized functions adapted to nonlinear 
calculations (and whose values are unbounded operators on a Hilbert space). 
These generalized functions were introduced 25 years ago, thanks to support of 
L. Nachbin [|71[8|| and L. Schwartz [9], but they suffered from the widespread belief 
that "there will never exist a theory of generalized functions adapted to nonlinear 
calculations." 

For the theory and applications of scalar valued nonlinear generalized func- 
tions, see the recent introductions in [|3l IH [HI, the survey [21 and the book [[TOl 
on their use in continuum mechanics, and the survey [|22| and the book [[l4| on 
their use in general relativity. The simplest definition of nonlinear generalized 
functions can be found in (22^ and references therein. The best and most recent 
textbook is [[Ml. 

One calculates with generalized functions as with C°° functions (although 
they seem basically very different, in particular when one attempts a physical 
interpretation): if ^7 is an open set in the generalized functions on are treated 
as if they were familiar functions that can be differentiated and multiplied freely, 
but there is a novelty. Classically if H denotes the Heaviside function it is natural 
to state = H. In the present context this is false: one has 7^ H. 

Let us explain this basic point in detail. In physics the Heaviside function 
H(x) represents a function whose values jump from to 1 in a tiny interval of 
length e around x = 0. Thus it is obvious that J (H^(s) — \l{x)).Lp{x) dx tends 
to when e ^ 0+ if (/? is a bounded function. But since H' is unbounded one has 
/ (H^(a;) — H(a;)).H'(a;) dx = 1/3 — 1/2 = —1/6, as obvious from elementary 
calculations. This shows that one is not allowed to state = H in a context where 
the function — H could be multiplied by a function taking infinite values such 
as the Dirac function 6 = H'. Therefore, in the ^-contextEone has to distinguish 
between "infinitesimally nonzero" functions such as — H and the genuine zero 
function, because "infinitesimally nonzero quantities" multiplied by "infinitely 
large quantities" can give significant nonzero results. For this reason we introduce 
a new notation: for two generalized functions Gi and G2 on an open set Q we write 
Gi ^ G2 ("Gi associated to G2") iff Vy? G C^{^) (i.e., infinitely differentiable 
with compact support in H) the number /^(G'l — G2) (x) (f{x) dx is "infinitesimal" 
(i.e., depends on a small parameter e > and tends to as e ^ 0+). 

'For brevity, we write "C/ -context" to refer to the nonlinear generalized functions context or 
setting. 
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The above was at the origin of widespread claims from mathematicians that 
"multiplication of generalized functions is impossible." Not only is it not useful 
and obviously impossible to keep classical equalities such as = H, but also 
physics displays instances in which several different Heaviside functions (asso- 
ciated to each other) are requested to model different physical variables: see the 
example of elastoplastic shock waves [5, fig. 6], [10, p. 106], O p. 185-189], 
[El [nil- From "concrete" examples in continuum mechanics P.l4l[5l[6l[T0ll2l[T8|| 
the basic mathematical idea became obvious: classical analysis was built at an 
epoch in which physics did not show "ambiguous products of distributions." The 
"ambiguities in products of distributions " that were found more recently in physics 
(in classical continuum mechanics [l3l|4l|5lll0l|T8l and references therein], in gen- 
eral relativity |[14l |22l and references therein]) appear to be due to an excessively 
strong mathematical idealization of physical objects which destroys ( or at least 
does not contain) the information needed to solve these problems. It is shown 
on physical examples in ['3l|4l[5l[Tp'| how nonlinear generalized functions contain 
more information. Indeed it was found that they were mathematically compelled 
to contain more information in order to solve the problem of multiplication of 
distributions [|7l[8l|9l, even before their use in physics began. 

To summarize one calculates with nonlinear generalized functions as in clas- 
sical calculations with infinitely differentiable functions (using =), but a basic 
novelty lies in that the classical equality of functions splits into = (the "true" 
or "strong" or "algebraic" equality) and the association ^ (also called "weak" 
equality). Concerning association the basic point is that Gi ~ G2 does not imply 
automatically G.Gi ~ G.G2 if G is another generalized function. But Gi ~ G2 
implies d/d^Gi ~ d/dr^G2. For example, ^ H but H^H' is not associated to 
HH': indeed f» H implies by differentiation 2HH' ^ H' and ^ H implies 
3H^H' H'. Therefore, association is some kind of weak equality not coherent 
with multiplication and strictly weaker than the (new) equality =. 

Vector valued generalized functions can easily be introduced, but the (presum- 
ably very difficult) problems stated in the abstract precisely come from the fact 
that we have to deal with generalized functions whose values are (densely defined) 
unbounded operators on the Fock space (i.e., a specific Hilbert space used in 
QFT). In continuum mechanics and general relativity one could use scalar valued 
generalized functions, what was far simpler and fortunately permitted to solve the 
problems with just nonlinear generalized functions. In QFT we are confronted 
with considerably more complicated mathematics. However, by making clear 
the mathematical background, nonlinear generalized functions permit a mathe- 
matical statement of the remaining difficulties (unbounded operators, computer 
calculations, and renormalization theory). 
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2 Difficulties in the physical interpretation of non- 
linear generalized functions 



Nonlinear generalized functions may have a (somewhat hidden) physical interpre- 
tation that is rather different from the one of classical functions. For instance they 
are not defined by pointvalues: two different generalized functions on an open set 
f2 in M may have the same pointvalues at all points of f2, [|l4l p. 37], see also IfTTll . 
The standpoint (of distribution theory) that the interacting field operators have to 
be smeared over a finite space-time region is not abandoned when we shall define 
them at a fixed point x = (x, t) of space-time by means of a nonlinear generalized 
function /. Indeed, because the "pointvalue" f(x) involves an "infinitesimal" 
neighborhood of x, the points have an "infinitesimal" but nonzero "extension," 
which is made to tend to in the definitions, presumably by ignorance of what it 
is in a physical context. 

Further, the physical interpretation of the word "infinitesimal" used above and 
the way to state physics with nonlinear generalized functions may depend on a 
scale (dimension, energy) of the phenomenon under consideration and may be 
based on an approximation at this scale: in [|3l IH [H the scale is the width 
of shock waves. Phenomena that have no classical analog occur in the ^-context 
even for solutions of simple partial differential equations (see Sec. 4] for the 
appearance of a void creating a "particle"). In short the concept of generalized 
functions compatible with nonlinear operations is extremely rich and subtle, with- 
out analogue in classical functions or distributions: there is an infinity of Heaviside 
step functions, of Dirac delta functions; generalized functions have an "internal 
structure," i.e., for each / each point x is in fact replaced by a portion of space-time 
(some kind of space-time "quantization?") that can be imagined as an approxima- 
tion of a Dirac delta function centered at x. Then x and / cannot practically be 
dissociated, and the pointvalue f{x) may even be related to the whole space-timejl 
thus introducing some very subtle kind of nonlocality, although at a rougher scale 
the locality expressed by the canonical commutation relations is fully satisfied... 

Therefore the underlying possible interpretations appear presently extremely 
numerous and one should avoid too hastly interpretations. In QFT the phenomena 
under consideration are much more mysterious than in classical mechanics, so one 
should be very cautious with interpretations El 

^Of course this "dependence vanishes very rapidly at infinity," remaining essentially, but not 
completely, a point dependence. 

^In contradistinction, the interpretation in 13] H] [5] |6l \W\ of the use of nonlinear generalized 
functions to model shock waves is rather clear because it reflects classical ideas of physics and it 
has been comforted by a number of calculations and experiments. 
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3 Aim of this paper 



This paper is focussed on the nonlinear calculations on distributions that constitute 
the basis of the Heisenberg-Pauli calculations to show that they can make sense 
mathematically. 

Wfe only pretend to give a mathematical sense ( very far from the " best one ") to 
the H-P calculations as a starting point for improvement. We do not try, consider- 
ing this would be premature, to interpret them physically. We stick mathematically 
to the H-P calculations in the simplest possible way: the interacting field operators 
will be usual functions of the time variable and generalized functions of the space 
variables because the time variable plays a privileged role in the (formal) H-P 
calculations^ Of course they should be presented in a next step as generalized 
functions on truly four-dimensional space-time. 

Problem: Modify the presentation below so that it becomes more suitable 
for the solution of the problems mentioned in the abstract, or for the physical 
interpretation. May be also consider more realistic physical models involving 
fermions. 



4 Prerequisites in QFT: The Fock space and the free- 
field operators 

4.1 Fock space 

The space of states called Fock space is the Hilbertian direct sum 

oo 

F = 0L|((M3)"), (4.1) 

n=0 

where, for n > 0, L|((R3)") is the Hilbert space of complex valued symmetric 
square integrable functions on (R'^)", with respect to the Lebesgue measure, and for 
n = 0, L|((]R^)") stands for the field of complex numbers. Thatis: an element ofF 

is an infinite sequence (/„)„, n = 0, oo, such that l/op+^lnridK/n)!!™)^ < oo, 
where || ||„ is the norm in L|((M^)") and /„ stands for the symmetric function 

^While the H-P procedure with its asymmetrical treatment of the time is not manifestly rela- 
tivistically invariant, it was shown by Heisenberg and Pauli that Lorentz-invariance is not destroyed 
by canonical quantization, what can also be verified by direct calculation, e.g., ||T6| . 
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(Ci)---)Cn) /n({i; Cn); 6 ^ Of coursc, from the definition of a Hilbcrtian 
direct sum, we have 

1/2 

ll(/n)l|F=(|/o|^ + E(ll(/n)y') ' (4.2) 

n=l 

+ 00 

((/n), (9n))F = /o ■ ^0 + X](/n, 9n)m (4-3) 

n=l 

where ( , is the scalar product in L|((M^)") and * is complex conjugation. In 
the sequel we shall use a dense subspace of F, i.e.: 

D = { {fn)n e F, such that = for n large enough } . (4.4) 

4.2 Creation and annihilation operators 

If e L|(M^) the creation operator a"'"(V') is given by the formula: 

a+(^) : (/n)n^ (o,/o ^(6),-, 

V^Sym(V^(e„)®/„-i),...), (4.5) 

where Sym is the symmetrization operator. If ^ G L|(M^) the annihilation 
operator a~ (t/j) is given by the formula: 

V^mn+l), /n+l(6, -,en+l)), ■■■)■ (4.6) 

The operators a+ {^) and a~ (ip) are defined at least on the dense subspace D of F, 
with values in D. They are not bounded operators on F because of the coefficients 
^/n and ^/n'+T. 

4.3 Scalar product 

In the quantum theory of a 'scalar' boson field the argument of the creation and 
annihilation operators is a Lorentz-invariant scalar function t^) e C where the 

variable {■Cj^^} = {^/i} G M^and/x = 1,..,4. The scalar product of two 1-particle 
states is then not given by the usual L|((M3)", C) formula (/i, ^i) = / d^^ f *^gi, 
and its obvious generalization to n-particle states, because the 3-volume element 
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d^^ is not Lorentz invariant. Instead, it is the relativistically invariant expression 

nasi p. 828], [[la p. 93] 

((/.(oii.(0)):^>///^_.^'«(/r|-f..). <^-'> 

which is a special case of 

«/.(oii.(0)):^.//X/E'-(/r||-f.). <^-«) 

where S is a space-like hypersurface, and fi a tensor index so that the contraction 
is indeed invariant. Taking for S the hyperplane = Cst leads to (|4.7I) . 



4.4 Free-field operators 

The free-field operator is defined by: 

(f>o{x,t)=a+{A+{^-x,t^-t))+a-{A+{^-x,^-t)), (4.9) 

where the variable {x, t} e corresponds to the space-time dependence of 4>q. 
The function 

A+(x,t) = — —e^vi{k-x-kh), (4.10) 

where fc" = \/k'^ + m? and A; G M^, is not in L|(M^). The mathematical object 
that makes sense is therefore: 

(Poi^,^) ■= 4>Q{x,t)(p{x) dx, (4.11) 

with Lp{x) a suitable test function. That is: (^q(x, t) is a distribution in the x- 
variable whose values 0o(v^' ^) densely defined linear unbounded operators on 
F (they map D into D). The conjugate free-field operator is: 

d 

7ro(a;,t) := —4)o{x,t), (4.12) 

which is again a similar distribution. Of course it would be better, in view of 
Lorentz invariance, to consider them as distributions in {x, t} but this is mathe- 
matically not necessary in order to stick to the H-P calculations. 
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5 Prerequisites in QFT: Sum-up of the Heisenberg- 
Pauli calculations in the simplest non-trivial case 



In this section the Heisenberg-Pauli calculations (see for example (24^, p. 21-22, 
292-332]) are summarized in the simplest non-trivial case, i.e., a self-interacting 
real scalar boson field. They are presented as purely formal calculations, i.e., they 
are not defined mathematically and are done by analogy with usual calculations 
on C°° functions. 

Let T,g eR. One sets ^ p. 206]: 

■^y^^^ i<M<3 

+ lm'{My,r)Y + j;^{My,r)f^'}d'y, (5.1) 

where 9^ = d/dx^. Of course this formula involves classically unjustified prod- 
ucts of distributions whose values are unbounded operators on a Hilbert space, and 
also an unjustified integration. Then one sets: 

(f){x,t,T) =exp(i(t-r)H(°)(r)).0o(a;,^)-exp(-z(t-r)H(°)(r)). (5.2) 

Again, the exponentials are not defined: they should be defined as operators (ex- 
pected to be unitary) on F, and there is also a problem of composition of operators 
because of the domain of (Pq{x,t). Then calculations mimicking calculations 
on C°° functions give (p{x, t, r) as the solution of the Cauchy problem called 
"interacting field equation: " 



—(f){x,t,T)= 22 -Q^(t>{x,t,T) -m^(t){x,t,T) - g[(t){x,t,T)) , (5.3) 



1<M<3 " ^ 

with the initial conditions 



d d 

0(x,r,r) = 0o(x,r), and — ^(x, r, r) = — (/)o(x, r). (5.4) 

It is a wave equation with nonlinear second member. Since the initial condition is a 
pair of irregular distributions the solution is not expected to be more regular than a 
distribution for which the nonlinear term does not make sense in distribution theory 
(with further a "big" problem due to the fact one is confronted with unbounded 
operators). Many such equations have been successfully studied in the ^-context, 
but in the scalar valued case, see the survey f,20.1 and references therein. 
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To solve (|5.3I) explicitly one introduces the operator: 

Sr{t) := exp(i(t - r)Ho) exp(-i(t - t)U.^^\t)), (5.5) 

where Hq is the "free-field Hamiltonian," i.e., Eq. (|5.1|) with (7 = 0. The operator 
Sr{t) is considered as unitary on F. Then, with, 

(f)j{x,t,T) = exp[i{t - r)Ho(r)).0o(a;,r).exp(-z(t - r)Ho(r)), (5.6) 

one has (from calculations): 

0(x,t,r) = (S.(t))-\0,(x,t).S.(t). (5.7) 

This formula suggests that the numerical results of the theory are the limits when 
T — oo, t +00 of the scalar products: 

lim |(($2||S.(t)<fi))F|', (5.8) 

where $i, $2 £ IP- One obtains (from calculations) that Sr{t) is solution of the 
differential equation: 

|:S.(t) = -^TTXT / {Mx,t,T)f^'d'x.SAt), (5.9) 
at iv + i Jxm^ 

Sr{r) = 1, (5.10) 

where 1 is the identity operator. This formula is the starting point for an attempt 
(called "perturbation theory") to calculate Sr{t) by developing it in powers of the 
"coupling constant" g, when g is small. 

Remark: These calculations form a rather intricate set of nonlinear calculations 
on distributions. Schwartz's impossibility result, see [|3]|, [fTOl p. 8], and [|f4l 
p. 6], has been interpreted as the proof that these calculations cannot make sense 
even if we forget that they deal with unbounded operators on a Hilbert space 
(which makes them considerably more intricate than calculations in the case of 
bounded operators, close to scalar calculations). But it appears now that Schwartz's 
impossibility result is circumvented (SecfH lIH IH [H [Ml ...]), so that we may give 
a rigorous mathematical sense to these calculations, although with severe needs of 
improvements due the fact that they deal with unbounded operators. 

6 Preparation 

If (y9 is a suitable function on we define its Fourier transform by: 

ifip)ik) = (27r)-3/2 / e-"'-''ip{x) dx. (6.1) 
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From now on Lp is such that its Fourier transform ?ip is C°° with compact support 
in and identical to 1 on a 0-neighborhood in M'^. 

In view of (|4.1 II) we first rewrite the free-field operator (14.91 ) so that it reads: 

0o(<^, e, a;, t) = a+ (A,(e -x,t^- t)) + a" (A,(e - x, - t)) , (6.2) 

where 

A,((^,e,x,t) = — — / — exp^(A;-x-A:°t).(F<^)(+eA;), (6.3) 

and e is a strictly positive real parameter that is arbitrarily close to 0. 

A further basic assumption will be that (y9 is real valued. Therefore 0o(v^' ^? ^) 
is a symmetric operator on F. We set: 

d 

7ro(¥?, e, X, t) = ^</>o(</') i^)- (6-4) 

00 (v^, e, X, t) and 7ro(v5, e, x, t) are densely defined symmetric operators (of domain 
D) on F. They map D into D and satisfy on D the commutation relations (notation 

[A, B]= AB- BA): 

[</)o(¥?, e, X, t), </)o(¥?, e, x', t)] = 0, (6.5) 
[7ro(v2, e, X, t), 7ro(v2, e, x', t)] = 0, (6.6) 
[0o(v2, e, X, t), 7ro(v5, e, x', t)] = i5{ip, e, x - x') 1, (6.7) 

where the nonlinear generalized function replacing the usual Dirac 5-function is 
given by 

<5(^,6,x-x') = ^(Fp)(^^), (6.8) 

where p is calculated from y9. The function p is C°° with compact support in R^, 
and is identical to 1 on a 0-neighborhood in M^. 

The above is equivalent to the classical description of the free fields considered 
as distributions, presented in a way convenient for the sequel. The formal descrip- 
tion of physicists comes about by taking (abusively) e = (or by formulating 
the free-field operators by means of the nonlinear theory of generalized functions, 
which for the sake of clarity is postponed to Sec.[IO]where it will be really needed). 
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7 Canonical Hamiltonian formalism: first step 



We set: 

+^"^^(</>o(<^.e,?/,r))^ + ^^q-y((/)o((^,e,y,r))^^^|.x(e?/) A, (7-1) 

with X a C°° real valued function on with compact support and identical to 1 in 
a 0-neighborhood|f| 

H(°) X, e, t) maps D into D and is symmetric. One can prove E p. 3 1 1-313] 
that it admits a self-adjoint extension denoted by H<*^>((/), £, e, r) on a domain 
denoted by containing D (do not forget that H<°^ and D<°> depend on 

V9, X, e, and r). Therefore, from the Hille-Yosida theory, {exp(itH<*'>)}tg]R is a 
strongly continuous group of unitary operators on F. We set (considering r fixed 
we do not note explicitly the r dependence of 0, similarly we leave the dependence 
on X implicit to simplify the notation): 

e, X, t) = exp(z(t - r)H<°>) .cj),{ip, e, x, r). exp{-t{t - r)H<°>) , (7.2) 

defined on 

©(t) := {exp{i{t - r)H<°>)©, (7.3) 

which is a dense subspace of F a priori depending on t (and also on cp, x, e, and 
r). We define 77(92, e, x, t) by formula (|7.2I) withTTo in place of 0o (we do not yet 
know that tt is the time derivative of (p). 

Simplification of the exponentials and the commutation relations (16.51 - 16.71) of 
the free-field operators give: 

[0((^,e,a;,t),</)((^,e,x',t)] = 0, (7.4) 
[n{(p, e, X, t), 7r(v9, e, x', t)] = 0, (7.5) 
[0(V9, e, X, t), 7r(v9, e, x', t)] = i5{ip, e, x - x') 1, (7.6) 

which are valid on D(t). 

^x(ey) is a damper introduced for the purpose of integrating over the whole of M.^. 



11 



8 The problem of domains 



From (17.21) . (f){(p, e, x,t) maps into D(t), and so does its expected time- 
derivative 



iH<°>. exp(i(t - r)H<°>).0o(<^, a;, r). exp(-z(t - r)H<°>) 

+ exp{i{t - r)H<°>).0o(¥', e, x, r). exp(-2(t - r)H<°>) .(-2)H 



In order to define the t-derivative of (f){(p, e, x, t) it would be desirable to have 
a domain independent of t, and this would be also desirable for the physical 
interpretation. (May be one should consider this 4> generalized function on 
space-time for this problem of domain.) 

Problem: Find an extension (/><o> of 4*0 such that cf) defined by (17.21) with 
^<o> would have a dense domain Di independent of t (and also if possible of 
x-i and r), or find another way to have a domain independent of t. (May be 
consider generalized functions in the variable (x, t) making Lorentz invariance 
evident.) 

In the absence of such a t-independent domain we content ourselves with 
a weaker definition of the t-derivative. We say that t ^ cf){t) from M to the 
set of linear operators from D(t) D(t) C F is differentiable with derivative 
t ^ r]{t) iff: for any map ^ : t ^ <l>(t) from M to F, with $(t) G B{t) and 
e D(t), Vt, then the map t ^ cf){t).^{t) from M to F is differentiable with 
derivative t i-^ ri{t).^{t) + cj)(t) (t) . With the above D(t) such a 77 is unique. 
Other definitions are possible to circumvent (at least provisionally) this domain 
problem. One obtains: 




(8.1) 



di 
d_ 

dt 



7r(v?,e,x, t)<l> = i[H 



7r(y?,e,x,t)]$, 



0((^,e,x,t)]$ 



(8.2) 



(8.3) 



for any $ in D(t). 
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9 Canonical Hamiltonian formalism: second step 



Set 

1<M<3 

+ (0(^, 6, y, t)Y + ^ 6, t))^-"'}. (9.1) 

The compositions of these operators, from D(t) into D(t), make sense. Moreover, 
since D(t) is independent of y, the 9^ derivatives make sense. Simplification of 
the intermediate exponentials yields: 

Hi^, e, y, t) = exp {tit - r)H<o>) .7^^ i^, e, y, r). exp(-^(t - r)H<o>) , 

(9.2) 

where H^"^ is defined like Ti. but with (/)g and ttq instead of and tt; Ti^^^ maps D 
into D. Again simplification of the exponentials gives: 

[ n{y^,e,y,t).xiey)d'y = -ti<'>. (9.3) 

Thus (l8:2tiO) give: 

d f 

—(t){if,e,x,t)<!> = i [n{ip,e,y,t),(f){ip,e,x,t)]x{ey)d^y<^, (9.4) 

— 7r((/3,e,a;,t) $ = W [n{ip, e,y,t),7z{ip, e, x,t)]x{ey) d^y ^, (9.5) 

if $ G B)(t). From (IT41 - I731) and the formula for H we get: 
d f 

—(j){(fi,e,x,t)^ = i S{ip,e,y - x).7r{(p,e,y,t)xiey) d^y ^, (9.6) 
^7r((/?,e,x,t) $ = 2 /" I ^ df,(f){ip,e,y,t).df,5{(p,e,y - x) 

Jj/GIR3 1<^<3 

+ "^Vlv', e, e, y - x) 

+ ^(0(y?,e,i/,t))^.5(y?,e,i/-x)|.x(ei/) d^y (9.7) 

if $ G D(t). At this point we are forced to use the nonlinear generalized functions 
that we had avoided up to now to provide a presentation involving only classical 
objects. These generalized functions provide a mathematical setting appropriate 
to take the limit e ^ 0+ so as to derive, in a suitable sense, the interacting field 
equation from (19.61 - 19.71) . They are also needed for the scattering operator. 
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10 The interacting field operators 



This section is only a brief summary. 

On D we define a "convenient kind of topological structure" f/l. p. 296^ and 
we define the algebra L(D) of linear operators from D into D that map any 
bounded set in D into (another) bounded set in D; by definition a subset M 
of L(D) is said to be bounded iff it is equibounded, i.e., \fB bounded set in 
D then {L{x);L E M,x G B} is again a bounded set in D; L(D) becomes a 
"bomological algebra" in the terminology of [7, p. 286] . Then we define an 
algebra Q(M.^, L(D)) of generalized functions whose values are linear operators 
from D into D [7, p. 287-295]. For fixed t, ^(R^ L(ro(t))) is defined by transport 
of structure through the exponentials that are in (17.21) . 4>{x, t) and 7r(x, t) are in 
Q (M'^, L(D(t))) , as generalized functions of the variable x G R^. 

Theorem ^ [Rigorous interacting-field equation] In the context of 
operator-valued generalized functions in Q one has, for all e 
C^, i.e., VH G C°° with compact support in the variable x, 

t/3xS(x)|-0(x,t)-7r(a;,t)} =0, (10.1) 

1<M<3 ^ 

-m^cP{x,t)-g{cj){x,t)f}=0, (10.2) 
with the free-field operators as initial values at time t = t. 



Translational and Lorentzian invariance: Let Ai be the Minkowski space and 
let a;, J G ^A. If T(d) is a translation and U (T{d)) its action on F, then translation 
invariance consists of checking the formula (l){x+d) = U {T{d)).cj){x).U (T(—d)). 
By following classical calculations, e.g., lfT6l . it holds in the present context taking 
into account adequate changes in r and x, even in the case of variable domains. At 
the end these changes disappear since r — > —oo and the function?/ — > x(ey) covers 
the whole space as e ^ 0+. The same applies to Lorentz rotations preserving the 
direction of time [|16I . 

''We only say that a set in D is bounded iff it is contained in a finite sum in the Hilbert sum 
defining F and each of its components is bounded; these bounded sets play the role of a topology 
although there is no topology but this is a natural extension of the concept of a normed space. 

'The proof of this theorem is given in Ref. lfT2l . 
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Problem (connected with the problem of domains): Improve the presentation so 
as to put in evidence mathematically the properties of the interacting field operators 
(hopefully they are "formally" correct since the above is a faithful transcription 
of the Heisenberg-Pauli calculations). In particular if = 3 formulate the 
properties listed in ['23', p. 96-102] (Wightman axioms), of course in an adequate 
mathematical way (due to the generalized functions setting and the results obtained 
so far), but by preserving their physical significance. 



11 The scattering operator 

We replace (|5.5I) by: 

S,(¥?,e,t) := exp(z(t - r)H<°>) exp(-i(t - r)H<°>), (11.1) 

where H^°^ £, e, r) is the self-adjoint extension of the interacting-field Hamil- 
tonian H'^°\ and Hq°^ that of the free-field Hamiltonian Hq which is equal to 
H<°> when g = 0. (To simplify the notation we leave implicit the dependence of 
Sr on X-) 

St-(v3, e, t) is a unitary operator on F. Then one has (after calculations): 

(f){(p,e,x,t) = {Sr{if,e,t)y\(l>j{(f,e,x,t).Sr{<f,e,t), (11.2) 
where ^/(v?, e, x, t) is the ^-version of (15.61) . i.e., 
(f>j{^, e, X, t) = exp(z(t - r)H<o>) .0o(<^, x, r). exp(-z(t - r)H<°>) . (1 1.3) 

The crucial step is of course the calculation of the t-derivative of St-((/9, e, t). 
We do this according to the definition introduced in Sec.[8l This leads to the 
rigorous differential equation 

d 

—Sr{(p,e,t) = ~iB.i{ip,e,t,T) S^(v3,e,t), (11.4) 
where the t-dependent interaction Hamiltonian is now 

Hj(<^,e,t,r) := ^ / d^x x{ex){(f)j{<^,e,x,t))^^^ , (11.5) 

in which the field operator 4>j{(p, e, x, t) is defined by (lll.3h . Note that (lll.5h is 
consistent because both H/ and cpj map © into ©. Therefore, in (II 1.41) . both Sr 
and dtSr map D(t) into D. 
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In combination with equations (II 1.21) and (|1 1.31) the differential equation (II 1.41) 
with initial condition S^((/9, e, r) = 1 solves the rigorous interacting field equation 
(I10.11110.2l) with the initial fields 4>Q{ip, e, x,t) and 7ro(v5, e, x, t) = dtcpoi^p, e, x,t) 
taken as those of a free field. 

In the limits r —oo,t +oo, which for now are formal limits, the operator 
St-(v5, e, t) = Sr{f, X; t) bccomcs the scattering operator, i.e., 

S(v5,X,e) := lim S_t(v?, X, e, (11-6) 

which still depends on the moUifier (p, the damper and the parameter e. The 
operator Sr{^, X; i) is unitary on F, and gives the transition probabilities (15.81) . 
That is, as "generalized real numbers," we have: 

V<l>i,$2GF |(($2||S.(<^,x,e,t)$i»p|', (11.7) 

but calculating their limits as r ^ — oo, t +oo, is non-trivial. 

Equation (111.41) corresponds formally to (15.91) . which leads to the (non-renorm- 
alized) perturbation series in the standard formalism. 

Problem: Is it possible to calculate S,- perturbatively, i.e., as a formal series of 
powers of g, from (II 1.41) ? Would something like the Ep stein- Glaser method work 
lfT3ll2T| ? Or else should one calculate globally starting from (|1 1.11) ? 

12 The problems of calculating the transition prob- 
abilities 

The transition probability that an initial state $i G D becomes after interaction a 
final state $2 G © is given by something intuitively looking like the limits: 

V{1^2)= lim |((<l'2||S.(t)$i))Fp. (12.1) 

T^ — 00 

Before taking the limits these expressions are "generalized real numbers" between 
and II $1 III ■ 11*^*2111 since the operator Sr{t) is a generalized unitary operator. 
The problem is to extract an "associated classical real number." The concept of 
association used till now in the context of nonlinear generalized functions consists 
of letting e — > 0+ and finding a limit independent of the particular function Lp in 
use, i.e., here such that F^p is C°° with compact support in and identical to 1 
on a 0-neighborhood in M!^. Since it seems doubtful that such a limit would exist 
a more general definition of association is needed. 
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Therefore, a natural probabilistic interpretation of the association (of a classical 
number to a generalized number) is introduced by an averaging process: if R is 
a generalized number of representative (v^, e) R{ip,e), intuitively bounded 
and oscillating as e — » 0+ (think for example at R{ip, e) = exp(cos(l/e)) which 
endlessly oscillates between 1/e and e as e ^ 0+) one may associate to ^ a mean 
value as e — s> 0+, provided it exists. For instance, one may associate to the above 
R the value 

lim / exp(cos—)dx= lim — / exp(cosx) dx, (12.2) 
.^0+ Jo e T^oo T Jo 

where the purpose of integrating over [0, 1] is just to capture the oscillations, as 
e 0+, in order to take their average. Other closely related averaging formulas, 
which give similar results, are possible. We choose (|12.2I) in the present context 
because it is the averaging formula considered in the theory of almost periodic 
functions [jTl. 

In Sr {(f, e, t) defined by (111.11) one has from the theory of nonlinear generalized 
functions to let e — 0+, and further in (|12.1I) one has to let r — > — cxo and t — > +oo, 
that is t — r — > +00. It therefore appears that at least the e-limit should likewise 
be dealt with by an averaging process as above. Since such an average always 
exists for all almost periodic functions one can ask: 

Problem: Is |($i, Sr,e(i^)'^*2) | an almost periodic function in the variable 1/e 
(or possibly in the variables t, r and 1 / e)? And if not, does a mean value exist in 
a natural sense as e — 0+ (or as e ^ 0+ and t — r — +oo)? 

From the theory of almost periodic functions, see the answer is clearly 
yes for far simpler cases such as | exp(iA/e)<l'2) where $i, $2 are in a finite 
dimensional Hilbert space and A is a symmetric operator. In our case we have to 
take into account the Hille-Yosida theory that is used to define exp (i(t — r)H<°>) . 

Problem: Try to calculate numerically such an average. 



13 Final comments 

We insist that the above should be considered only as an exercise on nonlinear 
generalized functions whose aim is to lay a groundwork for improvement — that is, 
as a convenient model for presenting a novel mathematical framework to improve 
the mathematical consistency of quantum field theory. 

In particular, our discussion has been restricted to the major steps going from 
the Heisenberg-Pauli formulation of the scalar field problem to a simplified defi- 
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nition of the scattering matrix. A comprehensive approach should discuss many 
pending problems related to the properties of the Green's functions, the defini- 
tion of asymptotic states, the derivation of closed-form and perturbative solutions, 
renormalization, etc. Moreover, the implications of nonlinear generalized func- 
tions should be related to the numerous efforts that have been made to give an 
axiomatic foundation to quantum field theory. 

An application of the nonlinear theory of generalized functions as done above 
should be complemented by deeper works that rely on other mathematical theories 
as requested in problems. After 1955 QFT has been influenced by the claim that 
"there will never exist a theory of generalized functions compatible with nonlinear 
operations," which implied that the Heisenberg-Pauli calculations, although they 
have always been the backbone of QFT [24J, would never make sense mathe- 
matically. Therefore the theory in |I71 [H |9J had to be first applied in classical 
continuum mechanics and general relativity (recent surveys |[3l|4l|5l|22l) in order 
to counterbalance this claim and be accepted by physicists as a possible tool in 
QFT, among other mathematical tools. 

The final conclusion is that the tool of nonlinear generalized functions — 
which is rather "soft" but basic to have a well defined mathematical background 
as starting point — appears to be convenient in that it opens doors by making 
mathematically understandable the overall framework of the Heisenberg-Pauli 
calculations. But, since QFT is very complicated, hard work using other needed 
mathematical tools such as the theory of unbounded operators on Hilbert spaces, 
and adequate computer calculations, are required before a final understanding can 
emerge. 

End note 

The detailed proofs of the calculations presented in this paper are given in Ref . 113 . 
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